Resonant Dimer Relaxation in Cold Atoms with a Large Scattering Length 
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Efimov physics refers to universal phenomena associated with a discrete scaling symmetry in the 
3-body problem with a large scattering length. The first experimental evidence for Efimov physics 
was the recent observation of a resonant peak in the 3-body recombination rate for 133 Cs atoms with 
large negative scattering length. There can also be resonant peaks in the atom-dimer relaxation rate 
for large positive scattering length. We calculate the atom-dimer relaxation rate as a function of 
temperature and show how measurements of the relaxation rate can be used to determine accurately 
the parameters that govern Efimov physics. 
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The Efimov effect is a remarkable phenomenon that 
can occur in a 3-body system of nonrelativistic particles 
with short-range interactions when the S-wave scattering 
length a is tuned to ±oo. Vitaly Efimov discovered in 
1970 that in this limit there can be infinitely many 3- 
body bound states with an accumulation point at the 
scattering threshold and a geometric spectrum []]]: each 
successive bound state is shallower than the previous one 
by a multiplicative factor. In the case of identical bosons 
with mass m, this factor is e 27T ^ s ° rs 515, where sq ~ 
1.00624. The binding energies of the Efimov states for 
a = ±oo can be expressed as 



E<j? ] = (e 27r / S0 ) n *- n h 2 n 2 Jm, 



(1) 



where tv « IS the binding wavenumber of the branch of 
Efimov states labeled by n*. The geometric spectrum is 
a signature of a discrete scaling symmetry with discrete 
scaling factor e*l SQ RJ 22.7. 

The Efimov effect is just one example of the universal 
phenomena characterized by a discrete scaling symme- 
try that govern the 3-body system in the scaling limit in 
which the range of the interaction is negligible compared 
to \a\ @, For a review of these phenomena, which 
we refer to as Efimov physics, see Ref. [J]. For simplic- 
ity, we focus on the case of identical bosons. We refer 
to the particles as atoms, their 2-body bound states as 
dimers, and their 3-body bound states as trimers. If the 
scattering length a is large and positive, there is a shal- 
low dimer with binding energy E D = h 2 /{ma 2 ). Efimov 
pointed out that there is an infinite sequence of negative 
values of a for which there is an Efimov trimer at the 
3-atom scattering threshold [2|]: a — (e 7r / s °)"a',, where 
a'„ w — 1.56/i^ 1 [5(. There is also an infinite sequence of 
positive values of a for which there is an Efimov trimer at 
the atom-dimer scattering threshold 0]: a — (e' r / s ") n o,, 
where a* w 0.0708^ 1 [5(. Another example of Efimov 
physics is an infinite sequence of positive values of a for 



which the 3-body recombination rate into the shallow 
dimer vanishes [a 0, @l • The universal aspects of Efimov 
physics are determined by two parameters: the scattering 
length a and the Efimov parameter k*. 

The alkali atoms used in most cold atom experi- 
ments have many deeply-bound diatomic molecules (deep 
dimers). Efimov physics is modified by the existence of 
the deep dimers, because Efimov trimers can decay into 
an atom and a deep dimer. Deep dimers also provide in- 
elastic scattering channels for three atoms or for an atom 
and a dimer. If there is an Efimov trimer near the scat- 
tering threshold, it can give a resonant enhancement of 
inelastic scattering processes. When a « (e^jV,, the 
resonant inelastic process is 3-atom recombination into 
a deep dimer 0,0- When a » (e*/ s °) n a*, the resonant 
inelastic process is dimer relaxation, in which the colli- 
sion of an atom and a shallow dimer produces an atom 
and a deep dimer [Toj | . If there are deep dimers, the 
universal phenomena associated with Efimov physics are 
determined by three parameters: a, k*, and a parameter 
77* that determines the widths of Efimov trimers [Io| . 

Finally, more than 35 years after Efimov's discovery, 
the first experimental evidence for Efimov physics has 
emerged. In a beautiful experiment with cold 133 Cs 
atoms, the Innsbruck group has observed a resonance in 
the 3-body recombination rate that can be attributed to 
an Efimov trimer near the 3-atom threshold [ll[. They 
used the magnetic field to control the scattering length 
of 133 Cs atoms in the lowest hyperfine state by exploit- 
ing a Feshbach resonance at —12 G. Since inelastic 2- 
body losses are energetically forbidden, the dominant 
loss mechanism is 3-body recombination. By varying 
the magnetic field from to 150 G, they were able to 
change the scattering length from —2500 oq through to 
+1600 ao, where ao is the Bohr radius. They observed a 
giant loss feature at a s» —850 ao- At the lowest temper- 
ature at which they measured the loss rate as a function 
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of a, which was lOnK, it can be fit rather well by the 
universal formula of Refs. 0, [H| with 77* = 0.06(1). 

The Innsbruck group also measured the 3-body recom- 
bination rate for positive values of a reached by increas- 
ing the magnetic field through the zero of the scattering 
length [ill ]. They observed a local minimum in the in- 
clastic loss rate at a ~ 210 clq. The minimum could be 
due to Efimov physics, although the value of a is not 
large compared to the van der Waals length scale for Cs 
atoms: (mCe/h 2 ) 1 ^ 4 « 200 Oo- Another complication 
was that the lowest temperature that was reached for 
positive a was 200 nK. Thermal effects can be taken into 
account if the 3-body recombination rate is known as a 
function of the collision energy. The universal predictions 
for nonzero energy have not yet been calculated. 

One way to confirm the discovery of Efimov physics 
would be to observe a second loss feature related to the 
first feature by the discrete scaling symmetry. In the case 
of identical bosons, the large size of the discrete scaling 
factor e 7 ^ w 22.7 makes this difficult. Since the 3- 
body recombination rate scales like a 4 , the loss rate at 
the second loss feature is larger than that at the first 
feature by a factor of 22. 7 4 « 2 x 10 5 . Another way is 
to observe the correlations between different aspects of 
Efimov physics that are predicted by universality. This 
requires the accurate determination of the parameters k* 
and 77* that govern Efimov physics. It may also require 
the calculation of the universal predictions at nonzero 
temperature to account properly for thermal effects. The 
effects of temperature on the resonant inelastic los s for 
negative a have been considered in Refs. 



the partial wave expansion 



For positive values of a, the Innsbruck group has also 
created mixtures of 133 Cs atoms and shallow dimers [lij ]. 
which has allowed them to study atom-dimer collisions. 
The most dramatic manifestation of Efimov physics for 
a > is the resonant enhancement of dimer relaxation 
for a near (e 7r/so ) n a,. The Innsbruck group has ob- 
served a resonance in inelastic atom-dimer collisions near 
a w 400 do [lij]. In this paper, we calculate the rate 
for inelastic atom-dimer collisions in a cold gas of atoms 
and dimers as a function of the scattering length and the 
temperature. By measuring the position and line shape 
of the resonance and using our results to take into ac- 
count the temperature, it should be possible to determine 
the Efimov parameters k* and 77* accurately. Once these 
parameters are determined, they can be used to make 
accurate predictions for other universal phenomena asso- 
ciated with Efimov physics. 

The collision energy for an atom and a shallow dimer is 
E = 3h 2 k 2 / (Am), where k is the wavenumber of the atom 
or dimer in the center-of-mass frame. The differential 
cross section for elastic atom-dimer scattering is 



d<7AD ,, n m ,2 
- 5r = l/A D (*,tf)| , 



(2) 



f AD (k,e) = j2 



2L+ 1 



^ fc cot 



ik 



P L (cos9). 



(3) 



The elastic cross section a ad (E) is obtained by integrat- 
ing over the solid angle of An. The leading terms in the 
expansion of k cot Sq D (fc) in powers of k define the atom- 
dimer scattering length a ad and the atom-dimer effective 
range r s ,AD' 



k cot 6£ D {k) 



'^/ a AD + \r s ^ AD k 2 + 



(4) 



The optical theorem relates the total cross section, which 
is the sum of the elastic and inelastic cross sections, to 
the 8^0 limit of the elastic scattering amplitude in 
Eq. ©: 



(total) 
7 AD 



An 



Im/ AD (M = 0). 



(5) 



If there are no deep dimers, the scattering is elastic up 
to the dimer-breakup threshold at E = Ed- 

We first consider the case in which there are no deep 
dimers, so the phase shifts S^ D (k) are real- valued below 
the dimer-breakup threshold. The phase shifts for L > 1 
are universal functions of ka only. The L = phase shift 
is sensitive to 3-body interactions at short distances and 
therefore depends also on a«*. Efimov used unitarity 
to derive powerful constraints on the dependence on a«* 
that he referred to as the radial law. He used the radial 
law to deduce an analytic expression for the atom-dimer 
scattering length a ad up to numerical constants that 
were first calculated by Simenog and Sinitchenko [III]. 
They have been calculated more accurately using an ef- 
fective field theory [B|, [l?], EH • This method was used in 
Ref. [5[ to calculate the atom-dimer effective range r s A D- 
These universal results are 



a-AD = (l-46 + 2.15cot[so ln(a/a*)]) a , 
r sAD = (1.13 + 0.73 cot[s ln(a/a») + 0.98])V 



(6) 



where fAD(k,9) is the scattering amplitude, which has 



Numerical calculations indicate that r s ^D vanishes 
within the numerical accuracy when a ~ (e 7r / s °) n 4.85 a* 
[j|. The above expression for r Si AD is compatible with 
Efimov's radial law and also satisfies this additional con- 
straint. 

We now consider the case in which there are deep 
dimers, so dimer relaxation provides inelastic atom-dimer 
scattering channels. The large binding energy of the deep 
dimer is released through the large kinetic energies of the 
recoiling atom and dimer. The event rate (3 for dimer re- 
laxation is defined so that the number of dimer relaxation 
events per time and per volume in a gas of atoms with 
number density ua and dimers with number density Ud 
is (3nAriD- If the atom and the deep dimer from the 
relaxation process have such large kinetic energies that 
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they escape from the system, the rate of decrease in the 
number densities is 



dt nA = dt nD 



-PUATlD 



(7) 



The event rate (3 can be expressed in terms of a statistical 
average of the inelastic atom-dimer cross section: 



k 



(inelastic) 
7 AD 



(E) 



(8) 



The inelastic cross section is the difference between the 
total cross section in Eq. (O and the elastic cross sec- 
tion obtained by integrating Eq. © over angles. If the 
temperature is large compared to the critical tempera- 
ture for Bose-Einstein condensation, the thermal aver- 
age indicated by the angular brackets in Eq. ([5]) can 
be carried out using Boltzmann statistics. In a gas 
consisting of atoms and dimers in thermal equilibrium 
(but not chemical equilibrium), the Boltzmann factor is 
exp(-3/i 2 fef/(4m/CBT)). 

In Ref. [lfj , we pointed out that in the scaling (or zero- 
range) limit, the cumulative effect of all the deep dimers 
on Efimov physics can be taken into account rigorously 
through one additional parameter 77* . The effects of Efi- 
mov physics can be dramatic only if 77* is much less than 
1. If the universal expression for a scattering amplitude 
for the case of no deep dimers is known as an analytic 
function of k* or a* , the corresponding result for a system 
with deep dimers can be obtained without any additional 
calculation simply by making the substitution 



In a* 



In a* — if]*/ sq . 



(9) 



Since the L = atom-dimer phase shift is the only one 
that depends on ft*, the L — channel is the only one 
that contributes to dimer relaxation at leading order. 
Making the substitution in Eq. ([9]) in Eqs. j6]), we obtain 
complex-valued expressions for the atom-dimer scatter- 
ing length and effective range: 



clad 



Ts.AD 



= (1.46 + 2.15 cot[s hi(a/a*) + w?*]) a, (10) 
= (1.13 + 0.73 cot[s ln(a/a*) + 0.98 + ^*]) 2 a . 



If the expansion for k cot 8q D (k) in Eq. ([4]) is truncated 
after the effective range term, the inelastic cross section 
reduces to 



(inelastic) 
7 AD 



(E) = 



4:Tr(-lmaAD - \\a A D\ 2 lYCLr St ADk 2 ) 



k\l + ia A Dk - \a A Dr s ,ADk 2 \ 



(11) 



Our truncation of the effective range expansion in Eq. ((4]) 
limits the validity of Eq. (Ill| to collision energies E <C 
Ed, where Ed = fi 2 / (ma 2 ). This expression could be ex- 
tended straightforwardly to the entire region E < Ed by 
using a parameterization of the atom-dimer phase shift 



below the dimer-breakup threshold given in Ref. The 
rate constant (3 for T <C Ed can be obtained by insert- 
ing the cross section in Eq. (fTT|) into Eq. ([8|). In the 
low-temperature limit, reduces to 67ra(— Im clad) /in, 
which can be expressed as [l(| 



= 



20.3sinh(2?7*) 



Ha 



sin [so ln(a/a*)] + sinh 77* m 



(12) 



If 77* is small, the maximum value of (3 occurs when a w 
[1 + sinh 2 r?*/(2s^)]a*, which is a « 1.0018 a* if 77* = 0.06. 
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FIG. 1: The dimer relaxation length m/3/h (in units of a*) as 
a function of a/a* for 77* = 0.06. The curves are for several 
values of the temperature: T = (solid line), T, /300 (dot- 
dashed line), T./30 (dashed line), where fcsT* = h 2 /(ma*). 

The resonant line shape for the dimer relaxation rate f3 
as a function of a is illustrated in Figs. Q] and [21 In Fig.QJ 
we show the dimer relaxation length m(3/Ti as a function 
of a/a* for rj* = 0.06 and several values of T. The tem- 
perature is given in units of T* = Ti 2 /(fc^ma 2 ), which is 
the dimer binding energy for a = a*. These results can 
be applied to any atomic system with large positive a by 
using the appropriate value of a* to convert to physical 
units. In Fig. [21 we show m(3/h as a function of a/a* 
for T = T*/30 and several values of 77*. As T increases, 
the height of the resonant peak decreases and its width 
increases, but the location of the peak remains very close 
to a* . If 77* is increased, the height of the resonance peak 
and its width both increase. The discrete scaling sym- 
metry implies that (3 has similar peaks at a ss (22.7)™a*, 
71 = 1, 2, . . . with maximum values of (3 that are larger by 
(22.7)™. The fact that the peak remains near a* implies 
that a* (and therefore the Efimov parameter k*) can be 
determined accurately by measurements at nonzero tem- 
perature without any need for thermal corrections. The 
value of 77* can be determined from the line shape of [3 as 
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a function of a, if the temperature is known to sufficient 
accuracy. 
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FIG. 2: The dimer relaxation length mf3/h (in units of a*) as 
a function of a/a* for T = T*/30, where fcsT* = fi 2 /(moJ). 
The curves are for rj, = 0.01 (solid line), 0.06 (dot-dashed 
line), and 0.2 (dashed line). 



hbach resonance? fs there any relation between the Efi- 
mov parameters associated with the large off-resonant 
scattering length and the Efimov parameters associated 
with the various Feshbach resonances? These questions 
can be answered by accurate measurements of the Efimov 
parameters. 

The Innsbruck experiment demonstrated that the res- 
onant inelastic losses associated with an Efimov trimer 
at the 3-atom threshold could be used to determine the 
Efimov parameters accurately for large negative a 



The Innsbruck group has also observed resonant inelas- 
tic losses associated with an Efimov trimer at the atom- 
dimer threshold [HI]. The results we have presented 
should make it possible to use those results to determine 
the Efimov parameters accurately for large positive a. 
These parameters could then be used to make universal 
predictions for other phenomena associated with Efimov 
physics, such as the locations of local minima and max- 
ima in the 3-body recombination rate. 

We thank R. Grimm for valuable discussions. This 
research was supported in part by the Department of 
Energy under grant DE-FG02-05ER15715. 



The resonance in inelastic atom-dimer scattering ob- 
served by the Innsbruck group occurs when the scatter- 
ing length is a w 400 ao [HI]. Since the van der Waals 
length scale for Cs atoms is (mCe/h 2 ) 1 ^ m 200 ao, this 
value of a is not deep into the universal region. If we 
ignore this problem, we can apply our universal results 
to 133 Cs atoms in the vicinity of the resonance by set- 
ting a* = 400 a . Using h 2 /(k B m) = 1.3032 Kag for Cs 
atoms, we find that the temperature scale set by a* is 
T* = 8145 nK. The temperatures in Fig. [1] correspond 
to 0, 27.2 nK, and 272 nK. The temperature in Fig. rj 
corresponds to 272 nK. The measurements of /3 by the 
Innsbru ck g roup were carried out at temperatures near 
250 nK (l5j. A reasonable guess for 77* is the value 0.06 
that was measured for large negative a on the other side 
of the zero of a. The value for j]* for this region of large 
positive a could be measured by fitting the line shape to 
our universal results, provided the temperature can be 
determined with sufficient accuracy. 

Universality requires that Efimov physics be described 
by the same parameters for large positive and negative a 
on opposite sides of a p oint where a — > ±00. The 133 Cs 
experiment of Ref. [11] studied Efimov physics for large 
positive and negative a on opposite sides of the zero that 
comes from the interplay between a broad Feshbach res- 
onance near magnetic field and the large off-resonant 
scattering length for 133 Cs. Is there any relation between 
the Efimov parameters for the two regions on opposite 
sides of a zero of the scattering length created by a Fes- 



[1 
[2. 

[3. 
[4 

[5 

[e; 
[7; 

p: 

[9 
[10 

in 

[12 

[13 
[14 

[IS 

[16 

[17 

[18 



V. Efimov, Phys. Lett. 33B, 563 (1970). 

V. N. Efimov, Sov. J. Nucl. Phys. 12, 589 (1971). 

V. N. Efimov, Sov. J. Nucl. Phys. 29, 546 (1979). 

E. Bra aten and H.-W. Hammer Phys. Rept. 428, 259 

(2006) [arXiv:cond-mat/0410417] . 

E. Bra aten and H.-W. Hammer, P hys. Rev. A 67, 042706 

(2003) [arXivxond-mat /020342T] . 

E. Nielsen and J.H. Macek, Phys. Rev. Lett. 83, 1566 
(1999). 

B. D. Esry, C.H. Greene, and J. P. Burke, Phys. Rev. Lett. 
83, 1751 (1999). 

P.F. Bedaque, E. Braaten, and H.-W. Hammer, Phys. 
Rev. Lett. 85, 908 (2000) [arXivicond mat 70002365] . 
E. Braaten and H.-W. Hammer, Phys. Rev. Lett. 87, 
160407 (2001) arXiv:cond-mat/010333l"] . 
E. Braa ten and H.-W. Hammer, P hys. Rev. A 70, 042706 

(2004) [arXivxond-mat /0303249] . 

T. Kraemer, M. Mark, P. Waldburger, J.G. Danzl, 

C. Chin, B. Engeser, A.D. Lange, K. Pilch, A. Jaakkola, 

H. -C. Nagerl, and R. Grimm, Nature 440, 315 (2006). 
JR. D'Incao, H. Suno, and B.D. Esry, Phys. Rev. Lett. 
93, 123201 (2004). 

S. Jonsell, Europhys. Lett. 76, 8 (2006). 
M.T. Yamashita, T. Frederico, and L. Tomio, to appear 
in Phys. Lett. A (2007) arXiv:cond-mat /0608542] . 
R. Grimm, talk at 18th International Few-Body Confer- 
ence, Santos, Brazil, August 2006. 

I. V. Simenog and A.I. Sitnichenko, Doklady Academy of 
Sciences of the Ukrainian SSR, Ser. A, 11, 74 (1981). 
P.F. Bedaque, H.-W. Hammer, and U. van Kolck, Phys. 
Rev. Lett. 82, 463 (1999) [arXiv:nucl-th/9809025] . 

P.F. Bedaque, H.-W. Hammer, and U. van Kolck, Nucl. 
Phys. A 646, 444 (1999) [arXiv:nucl-th/9811046| . 



